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Abstract 

In the spin boson model, the properties of the oscillator bath are fully characterized 
by the spectral density of oscillators J(ui). We study the case when this function is 
of Breit-Wigner shape and has a sharp peak at a frequency Q with width r <C O- 
We use a number of approaches such as the weak-coupling Bloch-Redfield equation, 
the non-interacting blip approximation (NIBA) and the flow-equation renormaliza- 
tion scheme. We show, that if f2 is much larger than the qubit energy scales, the 
dynamics corresponds to an Ohmic spin boson model with a strongly reduced tun- 
nel splitting. We also show that the direction of the scaling of the tunnel splitting 
changes sign when the bare splitting crosses $7. We find good agreement between 
our analytical approximations and numerical results. We illuminate how and why 
different approaches to the model account for these features and discuss the inter- 
pretation of this model in the context of an application to quantum computation 
and read-out. 

Key words: Spin boson model, quantum coputing, quantum measurement, cavity 

quantum electrodynamics 

PACS: 03.65.Yz, 03.67.Lx, 05.40.-a, 85.25.Cp 



1 Introduction 



The subject of open-system quantum mechanics and the physics of the bound- 
ary between classical and quantum physics has been of strong interest since the 
early days of quantum theory. A paradigmatic standard model for the study of 
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open quantum systems is the spin boson model [1,2]: A single two-state system 
coupled to a bath of harmonic oscillators described by the Hamiltonian 




(1) 



Here, H is a constant counter-term. The energy eigenvalues of the two state 
system alone are ±E/2 with E = \/e 2 + A 2 . The oscillator bath can model 
arbitrary Gaussian noise sources. It is fully characterized by a spectral function 
J{uj) which depends on the distribution of frequencies and couplings 



For a given physical system, e.g. a superconducting quantum bit coupled to a 
noisy electronic circuit, J(uj) can be obtained by analyzing either the effective 
friction [3] or noise [4] originating from the environment. It is useful to also 
introduce the semiclassical noise power S(lu) = J(uo) coth(cj/2T). 

Next to its long tradition in chemical physics, the physics of open quantum 
systems and in particular the spin boson model has gained recent practical 
importance in the field of quantum computation [5]. There, one is interested 
in obtaining long phase coherence times for the actual computation and long 
relaxation times for the readout. In a number of quantum computation real- 
izations, the researcher has the option to taylor or engineer at least part of 
the properties of the quantum system and the dissipative environment under 
study [6], e.g. in the case of superconducting qubits coupled to their con- 
trol and readout electronics [4]. In particular, environements with nontrivial 
internal dynamics, e.g. with resonances, can be realized and appear to be at- 
tractive [4,7,8]. Much is known about the physics of the spin boson model 
whose spectral density is a power law with an exponential cutoff [1,2]. Such 
spectral densities only contain the cutoff as an energy scale, which is typically 
assumed to be very high, leading to scale-free results. 

Much less is known about structured environments. We are interested in a 
generic realization of this physics described by a spectral density containing a 
Breit-Wigner resonance 
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in the underdamped case T f2. In that case, we can expand 

'(<*) = E ~i — — p n = n-2r*/n. (4) 



Moreover, we will be able to profit from the analytical continuation of results 
with a Drude spectral density, noting that eq. 3 can be written as 

J{uj) = % £ ™ . (5) 

v ; ait CT t^i uj 2 - {an + try v ; 



Note, that the shift of the resonance position from towards £7 can be ne- 
glected in the underdamped case, except close to the resonance. Hence, we will 
henceforth only emphasize this shift in those cases, when it actually affects 
the results. 

This type of spectral density is generically obtained by coupling the spin to 
a harmonic oscillator with eigenfrequency Q which in turn is damped with a 
linear friction coeffcient T/2. This friction is modeled quantum-mechanically 
by a bath of harmonic oscillators. Using a normal mode transformation, one 
can show that this is equivalent to our spin boson Hamiltonian with a struc- 
tured bath [9]. More details are given in section 5. This model is realized 
in various physical systems such as chemical reactions involving biomolecules 
[9], atoms in cavities [10] or superconducting qubits coupled to resonators 
[4,8,11,12,13,14]. It corresponds to the nonlinear dimer model of polaron physics 
[15]. The case of no dissipation with restriction to a rotating wave approxi- 
mation is known in quantum optics as the Jaynes-Cummings model [16]. Our 
notation corresponds to the one adopted in Ref. [7] and is slightly different to 
the one of Ref. [17]. 

The spin-boson model cannot be solved exactly and has been studied by a 
number of approaches. Some of them are largely numerical such as quantum 
Monte Carlo [18], real-time renormalization group [19], quasiadiabatic path in- 
tegrals [20] flow equation renormalization [21] and numerical renormalization 
group [22], others are mainly analytical such as the noninteracting blip ap- 
proximation (NIB A), a systematic weak damping approximation or exact Born 
approximation [23,24] or Bloch-Redfield [25,26,27]. A spectral density of the 
type eq. 3 poses a challenge to most of these approaches, since the spectral den- 
sity in units of the frequency, J(u)/u, can be either very small (off-resonance) 
or large (on- resonance) . In order to explore the phyiscs of this model and to 
obtain useful analytical information, these approximation schemes have to be 
applied within their range of validity and compared to numerical methods 
which are essentially nonperturbative in J(u)/u. Alternatively, one can treat 
the coupled TSS and oscillator system as multilevel quantum system and only 
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the friction to the oscillator as a bath [28]. 

The plan of this paper is to analyze this model using the weak-coupling Bloch- 
Redfield theory and the nonperturbative NIBA and to compare the results to 
a full numerical study obtained in the flow equation scheme. We will very 
briefly introduce these methods and compare the dynamics of the reduced 
density matrix [ characterized through the expecation value s z (t) = (a z ) (t) 
with, for definiteness, localized initial condition s z (0) = 1], effective reduced 
Hamiltonians, dephasing and coherence rates. Interpreteations of the results in 
terms of a superconducting quantum bit coupled to a resonant measurement 
circuit will be given. 



2 Bloch-Redfield 



The Bloch-Redfield-theory has originally been developped in the context of 
nuclear magnetic resonance [25]. It offers a systematic way to obtain a general- 
ized master equation within the weak coupling Born approximation between 
system and bath with J{uj)/uj as small parameter. It contains a subtle Markov 
approximation such that the resulting master equation is local in time; how- 
ever, the main bath correlations relevant within the Born approximation are 
kept and they do lead to time-dependent rates for a driven system [27,29]. 
Bloch-Redfield has been shown to be numerically equivalent to a full non- 
Markovian path integral technique for a rather generic choice of parameters 
[27]. Nevertheless, recent calculations at T = seem to indicate [24] that there 
may under certain circumstances be additional terms in the Born approxima- 
tion, that are neglected in the Bloch-Redfield approach. 

The natural starting point for the Bloch Redfield theory in the undriven case 
are the energy eigenstates of the spin-part of the Hamiltonian 1. In that "en- 
ergy basis", the Bloch-Redfield equation can be written as (see e.g. ref. [2]) 

p n m = —iuJ nm Pnm + ^ ^nmklPkl (6) 

kl 



where all indices take the values + and — corresponding to the ground and 
excited state and u nm = (E n — E m )/h. The Redfield tensor has the form 

J^nmkl u lm 1 nrrk ' °nk 1 Irrm 1 Imnk 1 Imnk V ' J 



where we have introduced 

T< imnk = (^z)i m ((rz)nkT(uj nk ) and Y\^ nk = (a z )i m (a z ) nk r*(-uji m ) (8) 
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where (cr z ) n k are matrix elements of a z in the energy basis, the * denotes 
complex conjugation. The basic building block of the rates in the Redfield 
tensor is the rate V which can be written as 

_^ oo oo 

T(5u) = — - / dte- i8ujt e-° +t [ dujJ(uj) [coscut coth(cu/2T) - i sincut] . (9) 
o o 



The resulting dynamics displays exponential decay and reads 

s z (t) = (e~™ + tanh (^) (1 - e'^)) + |f cos^Oe" 1 * (10) 



The quantities A cff and E cf[ = \Je 2 + A^ ff can be associated with a renormal- 
ized Hamiltonian 

H eS = \ (ea z + A eS a z ) . (11) 



The details of the shift of the tunnel splitting A i— > A c fj will be discussed 
below. 

The r r term in eq. 10 describes incoherent energy relaxation. It leads the 
system into thermal occupation of the renormalized Hamiltonian described 
below. The relaxation rate can be deduced from the Bloch-Redfield rates eqs. 
9, 7 and 8 



T r = R + R ++++ = (a z )_ + (a z )^ (T(E) + T(-E) + c.c.) 

= § 2 S(E). (12) 

This result is easily understood in terms of the Born-approximation: In order 
to relax, the system has to exchange the energy corresponding to the energy 
splitting E with the environment at once, using a single photon. 

The last term in eq. 10 describes quantum coherent oscillations analogous to 
Larmor precession of a spin [30]. These are the hallmark of (macroscopic) 
quantum coherence in the spin boson system. Their decay rate can hence be 
identified with the dephasing rate and can, using eqs. 9, 7 and 8, be written 

as 



r = -R e r 
r 



-+-+ - 



Re [2(a z 



r(o) + 



.(r(£) + r*(-£))] 

(13) 
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Note, that on very general grounds [30] we have 2T ( j > > T r . The extra factor of 
1/2 originates from the fact that there are in principle two dephasing channels 
corresponding to clockwise and counterclockwise Larmor precession. We are 
following here the standard NMR- motivated notation [30]; one could equiva- 
lently define 21^ as the true physical dephasing rate. The first term in eq. 13 
is proportional to the relaxation rate eq. 12, which reflects that a relaxation 
process certainly also randomizes the phase information. The additional term 
involves S(0), which in our case is oc T. This contribution originates from "fli- 
pless" dephasing processes which randomize the phases while keeping energy 
constant, i.e. transitions from a state into itself. 

The form of both rates eqs. 13 and 12 resembles the case of unstructured en- 
viroments [23], even though the spectral density eq. 4 has singularities close 
to the real axis. The high relaxation rate at E ~ f2 corresponds to resonant 
interaction between the qubit and the central environmental oscillator. When 
interpreting this result, one has to be aware, that the Born approximation 
involved is only valid for T r ^ <C E, which, bounding J(uj) < J(Q), means 
aVL 2 < T 2 . This a very rigorous constraint in the underdamped case, r <C Q, 
which we are considering. Also physically, we do not expect this result to be 
consistent up to strong couplings, because the relatively weakly damped big 
oscillator is a highly coherent quantum system which mostly reversibly ex- 
changes energy with the spin. However, since the Golden-Rule approximation 
in Bloch-Redfield only takes the long time limit, this reversible exchange can- 
not be seen in the Bloch-Redfield result. This can be understood from the 
order of limits prescribed by Bloch-Redfield and shown in the Appendix: The 
imaginary part of the energy is first sent to zero. Non-Markovian approxi- 
mation schemes [23,24] would at least take a self-consistent value and thus 
shift the S(E) in eq. 12 into KeS(E + iT r ). Such shifts can be important in 
particular if E ~ fi, when both predicted rates are very high. Results are 
summarized in figure 1. We clearly see the peaked behavior at resonance and 
notice that the influence of the self-consistent solution is rather small even 
at rather extreme parameters. On the other hand, the self-consistent solution 
predicts lower rates as compaed to the non-self-consistent one, similar to the 
predictions of flow-equation studies [17]. 

As mentioned above, the environment not only causes dephasing and relax- 
ation, it also renormalizes the tunnel splitting A (and with it the transition 
frequency), by dressing the two-state system with environmental degrees of 
freedom. This is similar to the physics of the Lamb shift or the Franck Con- 
don effect and leads, in the nonperturbative regime, to the dissipative quantum 
phase transition [1,31,32]. In our case, the transition frequency is renormal- 
ized according to E — * E — Imi? + _ + _. If we look at the imaginary part of 
the generic rate, eq. A. 2, F'(E) = -^r / dui J(u))V 2 ?L E 2 [coth(f3u/2)E — u] we 
observe a weight function V{u 2 — E 2 )~ l which changes sign at u ~ E. Thus we 
can expect an upward renormalization of E if most of the spectral weight of 
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Fig. 1. Relaxation and dephasing rates as a function of energy bias predicted from 
Bloch-Redfield theory. Parameters are A = 0.1S7, T = 0.01S7. For the low-damping 
plots we have chosen a = 10 -6 , T = W~ 2 Q, for the high-damping plots we 
have chosen a = 10~ 4 , T = 10 ft. The inset compares the self-consistent and 
non-self-consistent relaxation rate around the resonance for A = 0.1S7, T = 0, 
a = 1(T 2 and T = KT 2 fl 

J(uo) is above E (corresponding to E < f2) whereas E scales downward in the 
opposite case. Physically, this corresponds to level repulsion between the spin 
and the oscillators in the environment. The result also is consistent with usual 
second order perturbation theory for the energies. The sign changes happens 
at E ~ fi, the point where most of the spectral weight is concentrated, thus we 
expect a rather sharp structure of the splitting E eS (Q). Note, that this sign 
change is not predicted for the usual spin-boson in the scaling limit, which 
can be studied by the well-known adiabatic renormalization approach[l,31]. 
In that case, E e g is always reduced. This is consistent with our findings, be- 
cause in the scaling limit, the vast majority of the environmental oscillators 
have high frequency, much above the qubit splitting. 

From the structure of the dephasing rate eq. 13 we can conclude that the last 
term in eq. A. 3, which is even in energy, drops from the final result. Moreover, 
the remaining contribution to eq. A. 3 vanishes as E — > 0. If we finally go to 
low temperatures, we can replace p in eq. A. 3 by an appropriate logarithm 
and find for the shift of the transition frequency 



... A 2 a tQ 2 E x . 
oE = —- — 



E 2 2n T V E 2 - (aQ 2 + iY) 2 



log 



iE 



(14) 
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In the under damped limit we are working in, we can approximate the log- 
arithm as log|f2/i?| — icm/2 and split the result as E e g — E + 5E, 5E = 
SEfi + SE res . It contains a logarithmic contribution which resembles the scal- 
ing in the Ohmic case (with cutoff frequency Q), 



SEc 



2 

nE* 



J{E) log 



E 

n 



(15) 



This contribution changes sign from an upward shift at f2 < E to a downward 
shift at Q > E as is expected from the general arguments above. The logarith- 
mic divergence at low E can be indicated as a precursor of a dissipative phase 
transition. The other contribution takes into account the enormous spectral 
weight of the resonance, 

A 2 p2 _ Q2 

SE rcs = J(E)——. (16) 



This contribution is of the order a/Y. It will be shown below, that terms of 
this kind persist even in the absence of damping of the external oscillator. 
It, too, undergoes the expected sign change. It is linear at low E and hence 
does not contribute to a dissipative phase transition. It instead represents 
a substantial but finite renormalization. This is due to the fact, that for a 
dissipative phase transition, the environment has to get entangled with the 
spin down to arbitrarily small frequencies. These results are summarized in 
figure 2. As expected, we find in fig. 2, that the energy shift has a sharp 
structure around the sign change at E ~ Q. At this point, the spin becomes 
strongly entangled with the central oscillator, hence the concept of "qubit 
energy splitting" is of limited applicability. This observation is consistent with 
the usual dressed atom approach of cavity quantum electrodynamics. 



3 NIBA 



So far, we have restricted ourselves to the Born approximation, i.e. to the low- 
est order in J(uj)/uj and have otherwise kept the system general. We now turn 
to the noninteracting blip approximation (NIBA), which is nonperturbative in 
that parameter. It can be derived from evaluating the influence functional in a 
path-integral approach by assuming that the off-diagonal excursions ( "blips" ) 
of the density matrix contributing to the path of the two state system are 
uncorrelated [1,2]. It is thus justified when E <C fi, because then the bath is 
oscillating rapidly on the scale of the two-state system and the time-integrated 
bath correlation function quickly averages out, leading to weak damping on 
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Fig. 2. Renormalization of the energy splitting at T = taken at the degeneracy 
point e = 0, defined as positive if the splitting is decreased. Upper panel: Ohmic-likc 
logarithmic contribution from eq. 15; lower panel: Contribution of the environmental 
resonance from eq. 16 (for discussion see text). Low damping: T = 0.0117, a = 10 -4 , 
high damping: T = 0.1Q, a = 1CP 2 . 

longer time scales. Alternatively, the NIBA can be obtained by analyzing a 
the polaron-transfomed version of the spin-boson Hamiltonian. 

NIBA is known to work well under these conditions at the degeneracy point 
e = 0. At e ^ 0, the situation is more subtle. At e ^> A [1,2] the true dynamics 
is dominated by incoherent relaxation, which is again accurately predicted. 
This application of NIBA is closely related to the so-called P(E) theory of 
Coulomb blockade [2,33,34,35]. 

In this approach, the dynamics turns out to be governed by the Laplace trans- 
formed exponentiated correlation function 



P(A) = =- J e- xt e K ^dt 
o 



(17) 



where K(t) is the twice integrated bath correlation function from eq. 9 
1 °° J( ) 

K(t) = — [ dt-^p- ((cos cut - l)coth(w/2T) + ismcut) . (18) 
2n J ijj 1 
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At the degeneracy point, the dynamics of the system in Laplace space is readily 
found from 

oo 

s , W =/e- S «( i ) = XTW (19) 



where S(X) = (-P(A) + P*(A*))/2. Far from the degeneracy point, we find 
incoherent relaxation 



s z (t) = e 



-T r t 



1 



tanh I — — 
\2T 



+ tanh 



T r = 2ReP(ie + 0).(20) 



At T = 0, we can use eq. 4 to evaluate K(t) in closed form 



Kit) = ^5- V [e^- CT ' r )*Ei(-(^Q - aT)f) 

<7<t' 

-7 - log(-i((jQ - <rT)*)] (21) 

This is too compliated to allow a direct calculation of P(A) from eq. 17. At 
low energies, E <C f2 we can concentrate on the long time limit of eq. 21 and 
find, keeping only lowest order terms in Y/Q 

iWO = — jT" " 2 « Pog 1 0* | + 7 + wr/2] (22) 



where 7 is the Euler-Mascheroni-constant. This is a combination of a constant 
term of the order a /Y and a logarithmic term which resembles the findings in 
the Ohmic case [1,2]. This is similar to what we observed in our Bloch-Redfield 
result for the scaling in eqs. 15 and 16. From here, we find P(A) being 

A 2 / \ \ 2a 

P(A) = e -« n / r e - 2 T« e - <£ » r _r(l - 2a) - . (23) 



Off the degeneracy point, we can directly evaluate the relaxation rate from 
eq. 20 which reads 

„ A 2 n »e-^ /|e|\ 2a 

r '=ir' T iw fl • (24) 



This rate resembles to the Ohmic case [33,2] but is reduced by an extra ex- 
ponential prefactor exp(— aQ/Y), which again represents the contribution of 
the resonance and can be very small. Thus, we find te important result that 



10 



10 



-6 



0.001 



8=0. m 

8=1 0Q 




0.01 

a 



o.i 



Fig. 3. Relaxation rates calculated from the NIB A eq. 24, in the long time approx- 
imation using T = O.Olfi. Note, that for small a values the rate first grows with 
growing a, until the localization due to the resonance takes over and relaxation 
rates drop dramatically. 

by desigining small a but appreciable afl/T, the incoherent relaxation rate of 
the spin can be reduced to extremely small values. A physical interpretation 
of this finding will be given later on. The predictions of equation 24 are shown 
in fig. 3. At the degeneracy point, at e = 0, we find the Laplace transform of 
s z using eqs. 19 and 23. In analogy to the Ohmic case [2,1] we obtain for the 
backtransform that 



,(t)=E2-2a (-(A cff t) 2 - 2a ) 



(25) 



where E is the Mittag-Leffler function [2,1,36] and 



( e - an/r e- 27Q cos7rar(l - 2a)) 



l/(2-2«) 



(26) 



is the renormalized tunnel splitting. Note that this is only valid at A <C Q 
because we have taken the long time limit for K(t). Consequently, it always 
predicts a downward renormalization. As in the Ohmic case, the dynamics 
show a crossover from decaying oscillations at low a to incoherent decay at a > 
1/2 at e = 0. The renormalized tunneling frequency A c g shows a combination 
of the usual Ohmic scaling behavior governed by a, including a dissipative 
phase transition at a — 1, plus a very effective downscaling of e - a ^/$r(i-a) 
governed by a/F only, which also occurs for an undamped resonance and is 
not present for the Ohmic case. This again captures the contribution of the 
resonance and reflects the behavior we have observed in eqs. 16 and fig. 2 to 
lowest order in a/T. The dynamics is illustrated in figure 4. We can observe, 
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Fig. 4. Dynamical properties within the long time approximation of the NIBA at 
e = 0. Left panel: Effective tunnel matrix element eq. 26 for £1/ A = 10 for different 
values of the environment linewidth as a function of a. Lower T means more spectral 
weight and stronger downscaling of A c ff . Right panel: Expectation value of a z , eq. 
25 for different values of the damping parameters and A/17 = 0.1. 

that the time evolution of the spin is almost brought to a standstill, in the 
sense of absence of both oscillations and decay, already at modest coupling 
constants. 



4 Comparison to flow equation results 



So far, we have studied our system using traditional methods for open quan- 
tum systems. In order to complement this work, we want to compare the 
above results with previous work [17], in which the same setup was studied 
with the flow-equation renormalization method [21], which originates from 
strongly correlated electron systems and very well suited for treating prob- 
lems with several different energy scales. We will restrict ourselves to e = 0. 
This method typically can be used to calculate spin-spin correlation functions 
in equilibrium such as C(t) = (cr z (t)(T z (0)) ccl . A typical example is shown in 
fig. 5. The Fourier-transformed correlation function C{uo) is peaked at several 
frequencies. The resonance around A c g- corresponds to coherent oscillations, 
its width can be identified with the dephasing rate. There can also be a reso- 
nance around Q corresponding to oscillations of the oscillator leaving a trace 
on the qubit, but it hardly carries spectral weight. We have numerically solved 
the flow equations for small and moderate coupling strengths. More complete 
results are published elsewhere [17]. We see, that at A <C Q, A cS is rescaled 
downwards similar to the NIBA, but with quantitative differences The Bloch- 
Redfield result produces the correct slope at small a, see fig. 6. Around A = Q, 
the rescaling changes sign. Remarkably, Bloch-Redfield also predicts the slope 
above the sign change with good accuracy, see inset of fig. 6, although this set 
of date is taken very close to resonance. Please note, that in the inset fig. 6 it 
is important to keep Q in eq. 14. 
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Fig. 5. Correlation function as evaluated from the flow-equation method using 
fi/r = 0.06 and a = 0.15. 
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Fig. 6. Rescaling of the tunneling matrix element using different methods. Main 
plot: r/n = 0.02^, A/fi = 0.1; inset: T/Q = OMtt, A/fi = 1.1. 



5 Relation to quantum measurement and entanglement 



As already mentioned in the introduction, a straightforward way to implement 
this model with the spectral density eq. 3 is to couple the TSS to a single 
harmonic oscillator with resonance frequency Q, which is in turn damped by 
additional oscillators. This model has the Hamiltonian 
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H = —a? H a,r H 

2 z 2 2M 



2 +yfi 2 (X-g^) 2 



(27) 



The oscillator bath is characterized through an ohmic spectral density J(ou) = 

~2 

2m i uj- ^(- u} ~ Ui ) = ^Tu;, where, T/2 is the friction coefficient of the damped 
big oscillator. It was shown in [9], using a normal- mode decomposition, that 
this system is equivalent to the spin-boson Hamiltonian eq. 1 with spectral 
density eq. 3, where a = 2Mq 2 T/h. 

There are a number of realizations of such models. We would like to concen- 
trate on a realization in superconducting quantum circuits: A flux quantum 
bit coupled to the plasma resonance of a DC-SQUID. This setup has been 
thoroughly analyzed in Refs. [3,4]. It has been shown that the spectral density 
of the flux noise indeed leads to eq. 3 and how the circuit parameters relate to 
the parameters of that function. Moreover, it has been shown that the coupling 
parameter q actually can be tuned by the bias current through the SQUID. A 
representative circuit is shown in fig. 7. It is also shown there and discussed in 
Ref. [7], that a similar though less favorable realization can be found for charge 
quantum bits. We are mentioning this model, because it describes a detector 
of a quantum variable. Thus, we are going to interpret the findings of this 
paper in terms of quantum measurement theory. Other applications of reso- 
nantors coupled to superconducting qubits have been discussed in [12,13,14] 
As a key result, we have found above within Bloch-Redfield as well as within 
NIBA, that at e, A Q, the system dynamics can be interpreted as an Ohmic 
spin boson model with a strongly downscaled tunneling matrix element. This 
can be understood in terms of the following model, which was introduced and 
discussed already in Ref. [7]. We start from the undamped case, T = 0. the 
low-energy Hilbert space is spanned by |±) e ff = \^)\L/R) where |±) are the 
basis states of the qubit, o" z |±) = and \L/R) are coherent states of 

the harmonic oscillator centered around X = ±q, see Fig. 8. So in a general 
low-energy state = a\+) c s + b\— ) e s, \a\ 2 + \b\ 2 = 1, qubit and oscillator 
are entangled. In this low-energy basis, the Hamiltonian acquires the form 
of the renonrmalized spin part of the spin boson Hamiltonian eq. (11), with 
A c ff = A(L\R) = Ae~ v , where 77 = MQq 2 /h. This coincides with the result of 
eq. 26 in the limit of a — > but a/F = const. Under an appropriate choice of 
parameters, we can achieve i] > 1 and A e fr A. Following the notion of Ref. 
[37], the degree of entanglement is equal to 1 — e~ 2ri = 1 — |A e ff/A| 2 , i.e. we can 
interpret strong separation of the preferred states of the external oscillator and 
strong renormalization,i.e. A e fj/A 1 with strong entanglement. In terms of 
quantum measurement theory, the oscillator states are pointers onto the qubit 
states [38]. Chosing 77 ^> 1 corresponds to the condition of almost orthogonal 
pointer states in the environment, which has been identified as the condi- 



14 



Oscillator 



Meter 



R 




o 'b 





Flux-qubit 



charge-qubit 



Fig. 7. Underdamped read-out devices for superconducting flux (left) and charge 
(right) quantum bits involving a tunable Josephson junction representing a SQUID 
or a superconducting SET. 
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Fig. 8. Visualization of the ground state |0) and the coherent pointer-states \L) and 
\R) of the oscillator in the potential V(x) 
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tion for an ideal detector-dominated von- Neumann-measurement [38,39]. Such 
a measurement corresponds to the standard textbook quantum measurement: 
The preferred states into which superpositions are decohered are assumed to 
be Eigenstates of the measured observable regardless of the Hamiltonian of 
the qubit. In our case, eq. 27 describes coupling of the pointer degree of free- 
dom to a z and hence measurement thereof. Rescaling A c g asymptotically to 
zero means bringing the target states of the decoherence arbitrarily close to 
Eigenstates of a z , thus realizing the aforementioned textbook assumption. 

As it stands, the qubit just gets entangled with the pointers, but they are 
not read out. This can be done by coupling to the dissipative environment. 
As shown above, its influence corresponds to that of an Ohmic environment 
of strength a. Taking a < 1, this leads to dephasing and relaxation rates 
analogous to the Bloch Redfield results eqs. 12 and 13 



Note, that there may be nonexponential contributions to the dynamics as well 



It is important to notice that in the strongly entangled case, A c g <C A, the 
relaxation rate, which describes the thermalization of the system independent 
from the initial state, is strongly reduced, whereas the dephasing rate, which 
describes the projection of a superposition into a mixture of the eigenstates 
H e Q is hardly affected. This is a very favorable situation for a practical mea- 
surement: The information is quickly available, after = T^" 1 and is destroyed 
only after tr = T^ 1 . This is not only convenient for experimental implemen- 
tation but also guarantees high fidelity: The probability for reading out the 
correct result after the dephasing time is P = e~ T <t>l^ and thus close to unity. 
For completing the description of the measurement, one has to evaluate the 
resolution of the detector and the typical measurement times. This depends 
on details of the physical realization of interest and has been done in Ref. [7] 
for the superconducting setup. 



6 Summary and outlook 

We have studied the spin boson model with a structured bath using three dif- 
ferent approaches: Bloch-Redfield, NIBA, and flow equation renormalization. 
We have arrived at a number of common features: If the peak in the spectral 
density is at frequencies much above the environmental resonance, the system 
is equivalent to a renormalized Ohmic spin boson model. This has been inter- 
preted in terms of quantum measurement and the usefulness of this result for 





'eff 



(28) 



[24]. 
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modeling quantum detectors has been outlined. We have furthermore shown 
that the tunneling matrix element of the spin part is renormalized downward 
if its initial value A is below the environmental resonance Q and renormaized 
upward if it is above. We have compared this renormalization from all ap- 
proaches and shown that they are in reasonable agreement within the scope of 
their applicability. In particular, our analytical results from NIBA and Bloch- 
Redfield reliably approximate the numerical results from flow equations. 

This work would not have been possible without the extensive, illuminating 
and thorough work on the Spin Boson model by Prof. Uli Weiss. We would 
like to express our high regards and respect for him and his work and to 
congratulate him on his 60th birthday. 

We would like to thank M. Grifoni and S. Kehrein for useful discussions. Work 
supported by the ARO under contract Nr. P-43385-PH-QC and through SFB 
631. 



A Calculation of rates including poles of the spectral density 



We now want to outline how to calculate the rates eq. 9 We can interchange 
the order of integration and evaluate the time integral, which can be expanded 
into a delta function contribution and a Cauchy principal value. We can split 
T(E) into real and imaginary part, T'(E) and T"(E) and find 

F'(E) = ±J(E) [coth(/3£/2) - 1] (A.l) 

for the real part, which determines the decoherence and 

T ' {E) = 4^hJ duJ J{U)V ^& [coth(/3u;/2)£ - u] (A.2) 



for the imaginary part, which controls the frequency shifts. The latter can 
be calculated by extending the integration contour to the complete real axis, 
applying the residue theorem and resumming the resulting Matsubara series. 
We end up with 



„ a n 2 E ^ 



2tt 2iT ^ E 2 - (aVL + iYf 



p(T — icr) — Rep(iE) — ir 



E 



(A.3) 



where p(x) = ip{l + (3x/2n) + ip((3x/27i) involves the digamma function ip. 
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